Dynamical Evolution of a Cylindrical Shell with Rotational Pressure 
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We prepare a general framework for analyzing the dynamics of a cylindrical shell in the spacetime 
with cylindrical symmetry. Based on the framework, we investigate a particular model of a cylindri- 
cal shell-collapse with rotational pressure, accompanying the radiation of gravitational waves and 
massless particles. The model has been introduced previously but has been awaiting for proper 
analysis. Here the analysis is put forward: It is proved that, as far as the weak energy condition is 
satisfied outside the shell, the collapsing shell bounces back at some point irrespective of the initial 
conditions, and escapes from the singularity formation. 

The behavior after the bounce depends on the sign of pz, the shell pressure in the z-direction. 
When Pz > 0, the shell continues to expand without re-contraction. On the other hand, when 
Pz < (i.e. it has a tension), the behavior after the bounce can be more complicated depending 
on the details of the model. However, even in this case, the shell never reaches the zero-radius 
configuration. 

PACS numbers: 04.20.Dw;04.20.-q;04.70.Bw 



I. INTRODUCTION 

Investigations on how the spacetime structures would 
be during and after gravitational collapses are important 
topics in gravitational physics and spacetime physics. 
Nonlinear nature of gravity is making our understand- 
ing of the phenomena very difficult, both conceptually 
and technically. 

Steady development of numerical methods is contribut- 
ing to our deeper understanding on the phenomena. 
However, it does not mean the analytical investigations 
are becoming less meaningful. On the contrary, the more 
the numerical results we get, the more variety of ana- 
lytical investigations is needed to interpret them and to 
reach the whole understanding of the phenomena. 

In this respect, we have not yet obtained enough va- 
riety of mathematical models for gravitational collapses 
which can be studied, mainly by analytical methods. In- 
deed, the cases of just spherical collapses with various 
matter content is already complicated enough, and in- 
vestigations are still going on [1] . It can be inferred from 
this situation that constructing and analyzing other types 
of models may not be an easy task. 

As the next step in this direction, cylindrical collapses 
have been investigated considerably. This class of col- 
lapsing models are important because elongated matter 
distribution can drastically alter the destination of the 
gravitational collapse. Already in the Newtonian gravity, 
the Jeans scale plays a vital role and the elongated matter 
longer than its Jeans scale tends to split into smaller scale 
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objects. In the case of general relativity, there is the so- 
called "hoop conjecture" [2]; black holes are formed when 
and only when the matter of mass M is contained in a 
compact region whose circumference C in any direction 
always satisfies C < 2n ■ ^^r^- Therefore it is of great 
importance to study the fate of the collapsing elongated 
matter, such as cylindrical collapses. 

The numerical investigation of the collapse of the non- 
rotating dust spheroid by Shapiro and Teukolsky [3] sug- 
gested that the naked singularities could emerge when 
sufficiently elongated object collapses. This numeri- 
cal result motivated further analytical investigations on 
the cylindrical collapses. Here, in view of the cosmic- 
censorship hypothesis [4], investigations on more "nat- 
ural" conditions, rather than non-rotating cases, are of 
significance. In particular, the effect of rotation and the 
effect of matter pressure should be studied. Then, an 
analytical model of collapsing cylindrical shell made of 
counter-rotating dust particles was investigated by Apos- 
tolatos and Thorne [5]. (Here "counter-rotating parti- 
cles" indicates the situation in which the half of the par- 
ticles are rotating in the opposite direction to the other 
half of the particles.) It takes into account the rotation 
of particles through the rotational pressure, while the net 
angular-momentum of the shell is kept zero. By avoid- 
ing the frame-dragging effect in this manner, the analysis 
can be relatively simple. It turned out that even a small 
amount of rotation prevents the singularity formation as 
far as this particular model concerned. 

Another analytical model of cylindrical gravitational 
collapses was presented and briefly studied by Pereira 
and Wang [6, 7]. This model has interesting features in 
that not only has it taken into account both the rota- 
tion and the pressure effect, but it also tries to take into 
account the radiation of gravitational waves and mass- 
less particles from the collapsing region. Indeed, the 
numerical work suggests that the gravitational radiation 
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becomes important after the bounce of the cyhndrical 
shell [8]. In spite of several interesting features, however, 
this model has not yet been properly investigated so far. 

Since this model is what is investigated in this paper, 
let us briefly recall their work [6, 7]. Firstly they pre- 
pare general formulas in the coordinate system {t, r, z, (p} 
for matching arbitrary two cylindrical spacetimes. Here 
junction conditions [9] play the central role. Then they 
apply the formulas to look at a particular model. In this 
model, an interior flat spacetime and an exterior cylin- 
drical spacetime are matched together at r = i?o (t) ■ The 
geometry of the exterior spacetime is so chosen that there 
is out-going flux which can be interpreted as gravitational 
waves and massless particles (see Sec. Ill for the details 
of this model). Then a dynamical equation for the shell 
is expHcitly derived (see Eq.(39) in Sec. IV). Since this 
non-linear equation is complicated enough, a very spe- 
cial class of solutions is searched for, and three types 
of special solutions are claimed to be found. They are 
the solutions of (i) a collapsing shell which bounces back 
at some finite radius, followed by an eternal expansion 
afterwards, (ii) an eternally expanding shell, and (iii) a 
collapsing shell whose radius reaches zero, resulting in a 
line- like singularity formation [6, 7]. 

The solution (iii) looks in particular interesting, since 
it reminds us of the result of numerical calculations by 
Shapiro and Teukolsky [3], though the spacetime struc- 
tures in these two cases are not the same. Considering 
several intc^rc^sting features of the model, it is quite re- 
grettable that their investigation is a preliminary one in 
nature with the limitation to the very special class of so- 
lutions, and moreover it turns out that the derivation of 
the special solution faces with crucial problems, which 
could totally alter the main results (see Sec. IV for the 
detailed account on this point). Thus it is meaningful to 
analyze this model in detail. 

With these situations in mind, the present paper aims 
mainly two goals. 

Firstly, we set up a general framework for analyzing 
cylindrical shell collapses in the same spirit as Pereira and 
Wang [6], but paying more attention to the mutual rela- 
tions of three kinds of natural frames of reference. They 
are (1) the ortho-normal frame {e/i}^=o, 1.2,3 associated 
with the cylindrical spacetime geometry, (2) the ortho- 
normal frame {e„, e^, 62, e^} associated with a shell S in 
question (here e„ is the normal unit 4-vector to the shell 
S, while e-r is the 4-vector tangent to the world-line of 
an observer located on the shell), and (3) the coordinate 
frame {dt, dr, dz, d^}. Computations based on the frame 
(1) can be carried out systematically, while the frame (2) 
is more suitable for interpreting the results. Therefore 
the bridge between them would be helpful. The frame 
(3) has been also considered since it is often convenient 
for actual computations. 

The second purpose of this paper is to analyze the 
above-mentioned model, based on the framework we pre- 
pare, and investigate the dynamical behavior of the col- 
lapsing shell. We assume that the weak energy condi- 



tion [10] is satisfied outside the shell. It means, roughly 
speaking, that the collapsing matter behaves as normal 
source of gravity. It shall be proved that, the collapsing 
shell bounces back at some point irrespective of the initial 
conditions, and escapes from the singularity formation. 
The behavior after the bounce depends on the sign of 
the shell pressure in the z-direction, p^. When Pz > 0, 
the shell continues accelerated expansion, and its velocity 
asymptotically tends to the light velocity. On the other 
hand, when pz < (i.e. tension), the behavior after the 
bounce could be more complicated depending on the de- 
tails of the model. In either case, the shell never reaches 
the zero-radius configuration, contrary to the result (the 
solution (iii) mentioned above) of the preliminary analy- 
sis in Refs. [6] and [7]. 

In Sec. II, the fundamental formulas on the cylindrical 
spacetime and the cylindrical shell in it are derived. We 
in particular pay attention to the extrinsic curvature of 
the shell, which plays the central role for studying the 
shell dynamics. The jimction conditions, which is an- 
other essential part for the shell dynamics, are also inves- 
tigated there. Based on the framework so prepared, then, 
we investigate in detail a particular model of a collaps- 
ing cylinder in Sec. III. In particular, we prove a theorem 
that the shell never reaches the singularity in this model 
as far as the weak energy condition is satisfied outside 
the shell. Section IV is devoted to the detailed compar- 
ison between the present analysis and the one reported 
in Refs. [6] and [7], pinpointing why the latter had been 
led to the virtual singular solution. Section V is devoted 
for summary. 

II. GEOMETRY FOR A MOVING 
CYLINDRICAL SHELL IN A CYLINDRICAL 
SPACETIME 

A. The metric and frame-vectors 

1. The metric 

We consider a spacetime with cylindrical symmetry de- 
fined by a line-element, 

ds'^ = -T{t,rfdt'' + R{t,rfdr^ 

+Z{t,rfdz'^ + ^{t,rfd(j)^ . (1) 

(By a suitable coordinate transformation (t,r) (i;^); 
one can actually assume T = R, but we here keep the 
form Eq.(l).) 

Looking at Eq.(l), the theory of differential forms is 
efficiently applied [12], by introducing natural 1-forms, 

9" : = T{t,r) dt, 9^ := R{t,r) dr, 9^ := Z{t,r) dz, 
9^ : = $(i,r) d<t> . (2) 

The local ortho-normal frame {ca] a=o,i,2,z is then de- 
fined as the dual to the 1-forms {^"^}a=o,i,2,3- 
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Fundamental geometrical quantities, such as curva- 
tures, w.r.t. (with respect to) the frame {6^)^=0, 1,2,3 
are derived and collected in Appendix A for the conve- 
nience of later use and future applications. 

Now let S be a timelike 3-surface with cylindrical spa- 
tial symmetry embedded in a cylindrical spacetime de- 
scribed by the metric Eq.(l). Due to its symmetry, the 
surface E is characterized by a single equation r = p{t). 
Let r be a proper time of an observer located on the sur- 
face with dz = dO = 0. Then, the metric Eq.(l) yields a 
relation between dr and dt. 



dr' 



- (1 - |J p') de 



(3) 



We can use (t, z, (p) as a coordinate system on S, which 
defines ortho-normal vectors tangent to S, 



and 



(5) 



Here " ° " indicates the derivative w.r.t. r and X := 
dt/dr. 

A normal unit vector n of the surface S, which may 
also be written as e„ or 9„, and its dual one- form are 
given by, respectively, 



e„ = n = dn = {-P, -X, 0, 



= XTR{~ P /X, 1, 0, 0)(^trz 



{trzcp) 



(6) 



(The sufHx (trzcj)) indicates the frame employed.) We 
note that the set of vectors {en, e.r,e.z, e^} forms the local 
ortho-normal frame located on S. 



4. The symbols F, F' and F ° denote the partial 
derivatives of a fimction F w.r.t. t, r and r, re- 
spectively. In the same manner, dnF indicate the 
directional derivative of F w.r.t. the vector h. 



3. Frequently used formulas 

We here summarize essential formulas frequently used 
below. We note the relations, 



dt _ 
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o dp . dt . 
, P:^-^=p — =Xp. 



Therefore it follows 



(7) 



(8) 



For any function F(t, r) defined in the neighborhood of 
S, it follows 

F = XF+ PF' = X{F + pF') , 

dnF = T-^Rp +XTR-^F' 

= XT-^R{pF + T'^R-'^F') . (9) 

Conversely, we also have 



F = XT^ F-TRp dnF , 
F' = -R^ P F+XTRdnF 



(10) 



Each derivative in Eqs.(9) and (10) is understood to be 
estimated on S. 



Extrinsic curvature of the surface S 



2. Summary of notations 



1. In the ortho-normal frame { }a=n,T,z,4 



At this stage, let us summarize notations employed 
throughout this paper. 

1. Capital Latin indices A. B, C, • • • arc reserved 
for the indices w.r.t. the local ortho-normal frame 
{eA}A=o,i,2,3, and runs from to 3. Einstein's 
summation convention is adopted throughout the 
paper. 

2. The first few Greek letters a, (3, 7, • • • are re- 
served for the indices w.r.t. the ortho-normal frame 
{in, Br, ez, e^}, and take the values n, r, z or ^; The 
later Greek letters - ■ ■ are used for the indices 
w.r.t. the coordinate bases {dt, dr, dz, d^}. 

3. Since several frames appear in the paper, we indi- 
cate the frame explicitly when necessary. For in- 
stance, {p,q,r,s)(^trz4>) ■■=pdt + qdr + rdz + sd^{oi 
a vector and (p, q, r, s)(tj.z(t>) •= pdt+qdr+rdz+sd^ 
for a 1-form. 



We now compute the extrinsic curvature of S, which 
later plays the essential role in analyzing the shell dy- 
namics. The extrinsic curvature is handled mainly in the 
oiiho-iwrrnal frame {e-r, e^, e^} in this paper. To remind 
us of it, we write i^e^e^ etc. rather than Ktt etc. 

There is a standard formula which relates K^^e^ with 
the Christoffel symbol T and the normal vector h, 



K, 



eae/3 



■ a/3 

= —E^ cxE^ c^3,B ~ E^aE^'sd^E^ 13 ,(11) 

where in the last line, Eq.(Bll) in Appendix B has been 
used to express the Christoffel symbol T'^^^ in terms of 

since the latter, given in Eq.(A3), is much easier to 
compute. 

Let us first compute Ke^e^, 

Ke^e-r = —E^rE^rE^C^^B ~ E'^rE" nd^E^ ^ ■ (12) 
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The final result is 

i^e.e. = -^P-J^P(lni?Y + i(T^X^-l)a„(lni?^) 

+^^(X2t2 - 1) dn HT^/R^) - P {X^T^ - 1) (ln(rVi?2))° . (13) 



We note that Kc^c^ contains the second time-derivative 
of p as is seen in Eq.(13), so that K^,^,,^ plays the most 
important role in the shell dynamics. In view of its es- 
sential significance, the detailed derivation of Eq.(13) is 
shown in Appendix C. It is by far easier to compute 
Kf,^^^ and i^e^e^- Combining Eq.(ll) with Eqs.(B6) and 
(BIO), we get 

i^e^e, = a„lnZ . (14) 
i^e,e, = a„ln$ . (15) 



Here the r-derivative in the last term on the R.H.S. 
(right-hand side) has been purposefully left untouched 
to make the formula look simpler; the term will vanish 
whenever T = ii, as is the case for the model we analyze 
later. 



C. Junction Conditions 

Discontinuity in the extrinsic curvature across a 3- 
surface implies that the energy-momentum is accumu- 
lated on the shell; this is analogous to the relation of the 
discontinuity in the electric field across a shell with the 
electric charge there. 

Geometrical discontinuity across a 3-surface E is de- 
scribed by a discontinuity-tensor induced on S defined 
as 

KSo,0={K^0\-h^0{K\ , (17) 

where [Q] := (5+ — Q-, the difference of Q across the sur- 
face; K := SttG is the Einstein's gravitational constant. 

Then the discontinuity is governed by the junction con- 
ditions [9], 

(I) ds^_^^=dsl^^ (II) D^Sj = [T^^E''„n''] . 

(18) 

Here T^j^ is the energy-momentum tensor in the neigh- 
borhood of S; and Ef^a are given in Eqs. (6) and 
(B6), respectively, if one refers to the {t, r, z, </))-frame. 



I 

2. In the coordinate frame { }n=t,r,z,<i> 



It is reasonable to compute every quantity in the ortho- 
normal frame {ea}a=n,T,z,4> as has been done just above. 
However, in order to compare the present results with 
the results in the preceding work Ref. [6, 7] (see Sec. IV), 
it is also needed to represent main quantities in terms of 
the coordinate frame {9^}/i=t,r,2,0- 



(16) 

I 

Let us now prepare more explicit expression for 
Eq.(18). The metric induced on the surface S is 

= -^^ r O O , 

where 0'^ = dr, 9^ = Zdz, and d"^ = <^d(j) form a set of 
ortho-normal 1-forms on S dual to the bases {ea}a=T,z,<t>- 
Thus the discontinuity-tensor in this frame becomes 

rv 

5e,e, = ([i^e.e.] - [i^e.ej) =: , (19) 

where e, Pz and are interpreted as the energy density 
of the shell, the pressure in the ^-direction and pressure 
in the (?i-direction, respectively. It is easy to get 

Tli,^ = {InZr , r^;,^ = (ln$)° , T^^,^ = (hiZ)° , 
Te^e. = (ln$)° , Others = . 

(We note {ere^e^} forms the non-coordinate bases and, 
in particular, r!% = T^*- = 0.) Thus 

D0Sl^^ = {- e-{e + p,){\nZy-{e+P4,){ln^r, 0, 0) 

(20) 



I 

After some computations (see Appendix C for more details), we finally reach the result, 

i^e.e. = -X^TR p + {X^T^ - 1)^ (^(InT^)' + |J p (Ini^^) ^ + ^X^TR p {HT^/R^)y 
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in the (^^^;2:^^)-coordinate. Then (II) in Eq.(18) reduces 
to 

e +(e + p,)(ln Zf + (e + p^)(ln $)° = - E\ , 

(21) 

with n'^ and E^a being given in Eq.(6) and Eq.(B6), 
respectively. 
Noting that 



s.. = er o r + 33^61^ o r + p^6i^ o r 

= edr dr + p^;Z{Tfdz rf^ + p^^{T)'^d(l) d(j) 

We also get 

Sab = diag{€, Z(t)^p^ , $(t)^P0)(^,^,^) . 

It is straightforward to get 



1 



1 



= ^{zy , r^, = -($^)° , r^, = ri^ = (inz)° , 

rt^ = rt = (ln^r , others = , 
and 

DtS''a = {-e-{e + p,){\nZr-{e + p^){\n'i>y, 0, 0) 

in the (t, z, (/))-coordinate, which is the same expression 
appears as in Eq.(20). Therefore, (II) in Eq.(18) reduces 
to Eq.(21) in the (r, z, (?i)-fraine also. 



III. ANALYSIS OF A COLLAPSING SHELL 
MODEL WITH ROTATIONAL PRESSURE AND 
GRAVITATIONAL RADIATION 

Having prepared all the necessary formulas in the 
preceding sections, let us now reanalyze the cylindri- 
cal collapsing-shell model considered by Pereira and 
Wang [6]. The spacetime to be analyzed is constructed 
by matching two spacctimcs A4± described by the met- 
rics ds± at the timelike surface S (the suffix — is for the 
interior geometry and + is for the exterior geometry): 
The ds± are given by 

ds'i = -dtJ +drJ + dzJ' +rJd(f)J, (22) 
ds\ = e2T«)(-dt+2 + rfr+2)+d^+2+r+2d(/>+2(23) 

where 7 = 7(^), a function of ^ := t+ — r+ only: The 
surface E is assumed to be described by r± = p±{t). 
The metric ds+ has the form of ds in Eq.(l) with 

T{t+,r+) = R{t+,r+)=e^^^\ Z{t+,r+) = l, 
^{t+,r+) = r+ . 

The C-energy (cylindrical energy) [5, 11] of this space- 
time is given by 



(24) 



SO that we assume 7 > from now on. 

The interior spacetime is obviously flat, while the ex- 
terior geometry is given by, with the help of Eq.(A8) and 
noting that 7 = — 7', 7 = 7", 

R-oo = R-ii = ~e '^^ , 



r+ 

7 

— ( 

r+ 



W-oi — -"^lo — — e — e 



others = 



Therefore R = and Goo = Gn = — Goi = — Gio = 
or 

G., = ^€-'^^0'^ (8)9° -e° (^e^ -e^ (g)e° + (^e^) 

r+ 
7' 

= — {dt+ dt+ — dt+ dr+ — dr+ dt+ 
r+ 



+dr+ rfr-f 



Thus [6] 



'ab 



kakb 



(25) 



where A:^ = (a, —a, 0, 0)(^trzct>') is a null- vector. 

Combined with the Einstein equation, Eq.(25) implies 
7' > when the weak energy condition for matter is 
imposed. Since we are interested in the behavior of the 
collapsing shell in more or less "natural" conditions in 
view of the cosmic censorship conjecture [4] , it is reason- 
able to assume 7 > and 7' > from now on. 

It may be appropriate to make some remarks on the 
energy condition. Without any condition imposed on 
matter, any strange spacetime would become possible. 
Indeed, one can always claim that arbitrary spacetime 
is the solution of the Einstein equation if the matter 
with energy- momentum tensor Tab '■= '^^ab is prepared. 
Therefore it is essential to assume some kind of energy 
condition to make any meaningful arguments [10]. The 
weak energy condition is one of such appropriate condi- 
tions widely accepted, so we assume the condition here 
to make the definite arguments below. Of course, there 
is still a room of consicicring the violation of the energy 
condition at the shorter scale, due to the Casimir effect 
for instance. Though such a possibility is certainly in- 
teresting to study, we do not consider it in this paper 
for definiteness and only investigate the purely classical 
situations. 

We now impose the junction conditions Eq.(18) for two 
geometries A4±: 



Junction Condition (I) 



We impose 
dsl 



r-=P-(t-) 



ds\ 



(26) 
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When dz± = d(f>± = 0, it implies 



de_-dri = {i-pi)de_ 

e'''^+{dt\-dr\) = e''^+{\- pl)dt\ 

where 7+ := 7(t+, By plugging the relation 
pi = (3!^) P+ iiito the above relation and by solving 



Let us find the expressions for i^'e^Cb (Eqs.(14), (15) 
and (16)) in this model, and then derive the formulas for 
e, andp^ (Eq.(19)). 

First of all, the factor X becomes 



X 



the latter w.r.t. we get an important formula 



dt+ 
dtl 



{(l-e2^+)p2^e2^+}-'/'=:A 



(27) 



implying that 



Taking the f-derivative of the above relation, one finds 



dH 



1 



+ -^{p+(l-e2^+)p+-e2'^+7;(l-p^)} • 

(28) 



de 



Then, noting the relation X_ = AX+, Eq.(16) with 
E? = e^T yields 



Noting the relation. 



d 



1 



A = -p+p+-— (p+p++7;(l-p^)) , (29) 



dt+ A 
one also finds 



P-(M = A-^(A-V.) 



^ = -X+'e''r+p++j'^X+'e''^^pl{l-p+) , 
_ = -XJp. 

where Eq.(30) has been used in the last line. Noting also 



e27+ 



{p++^'+{1-pI)P+} . (30) 



K. 
K, 



= (9„ In r = — 
r 



Junction Condition (II) 



one finds 

Equation (19) then yields 
Ke ~" 



-7+ 



A(l-p^)3/2 



{(A - l)p+ - 7;p+(l - p+)[(A + l)p+ + 1]} 



= 



A - 1 



p+^l-p\ 



A - 1 











A(l-p2_)3/2 




{(A- 



(A - l)p+ - 7;p+(l - p+)[(A + l)p+ + 1] - A(A - 1) 
{(A - l)p+ - 7>+(l - p+)[(A + l)p+ + 1]} , 



P+ 



(31) 

(32) 

(33) 



with A = [pX + e^-^+il- p\)\ 



1/2 



Noting the relation A^ — 1 = (1 — p^){e^'^ — 1), we see When 7 = 0, A = 1 and e = pz 

(34) 



0, as it should 



> 

7+ = 

< 



> 

A = 1 

< 



> 
e = 

< 



7 



be. 



Setting 



Eq.(32) yields the dynamical equation of the shell, 



Po := Kp^ 



p+ = A 



■(l-p+2)=^/V+Po 



(35) 



To make the statements concise, we indicate the first, the 
second and the third terms on the R.H.S. of Eq.(35) by 
the symbols, [[!]], [[2]] and [[3]], respectively. 

The term [[3]] indicates that when the shell has a pos- 
itive pressure in the z-direction (Pq > 0), it causes more 
acceleration in the expanding direction, while the tension 
of the shell (Pq < 0) works to reduce the acceleration. 
However the term [[3]] does not change the whole dy- 
namics so much: The term could become important only 
when A ~ 1, i.e. when ^ ±1. Setting /?+ = ±(1 — (5) 



(i5 > 0), the term behaves as 



that the con- 



tribution of the term is not significant. Indeed, we shall 
show below a theorem that under the dynamical equation 
Eq.(35), the shell never reaches the singular configuration 
p+ = irrespective of the initial conditions nor the value 
of Po. 

First of all, we note the following obvious fact: 
Lemma 1 



Consider the case Pq > 0. Once the shell satisfies 
the condition > 0, it expands forever afterwards. 

Proof: 

Noting that A > 1 and that the terms [[2]] 
and [[3]] are non-negative when p+ > 0, it 
1— 

follows that p+ > ^ > 0, so that the 



claim follows. 



□ 



In particular, when \p+\ <C 1, Eq.(35) becomes p+ ~ 
e'^+/p^, so that the effective potential behaves like 
Ceff ^ —e'^+ log p^ + const, in the region where 7+ 
does not change so drastically. Compared with the dust 
case [5], where Ce// ^ a/p"^ + b\np_^_ + const, (a, b 
arc constants), the tendency of expansion of the present 
model is obvious. 

Next, let us investigate the contracting shell. To avoid 
any confusion, let us set V := — p+, so that < ^ < 1. 
Then Eq.(35) reads 



P+ = A^— ^ + 7;^(1 + V)[{A + 1)V - 1] + ^(1 - V'f/'e-r+Po , 



with A = [V^ +e^-^+{l-V^)] 



P+ 
1/2 



(36) 



r 



The present model has been constructed by giving met- 
rics Eqs. (22) and (23), and there is no particular physi- 
cal image regarding the matter content forming the shell. 
However, one can infer that the term [[!]] comes from the 
angular momentum effect since it is a decreasing function 
of p+. On the other hand, the term [[2]] may be identi- 
fied with what corresponds to Newtonian gravity under 
suitable conditions, since it is proportional to 7^. 

Let us pay attention to the last factor in the term [[2]], 
f{V) := (A + l)V - 1. The function f{V) is obviously 
monotonic, increasing, continuous function of V , with 
/(O) = —1 and /(I) = 1, so that fiy) has only one zero 
in (0,1). 

When the weak energy condition is satisfied, and the 

shell is not moving so fast (V^ ^ 1), then f{V) < 0, so 
that the term [[2]] yields an attractive force, which is con- 



sistent with the attractive nature of Newtonian gravity. 
A curious feature arises, however, when V ^ 1. Then 
f{V) > 0, and the same term now yields a repulsive 
force. Due to this characteristic behavior of f{V), the 
shell bounces back at some point, and never collapses to 
singularities as we shall see soon. 

Before proving this statement, let us first get rough 
understanding on what happens when p+ becomes small. 
The term [[!]] is the repulsive force, produced by an infi- 
nite effective potential barrier near p+ = (angular mo- 
mentum effect). In case the shell would have ever reach 
/9+ = 0, it should have been this repulsive effect becomes 
negligible compared to other terms, i.e. when p'^ ap- 
proaches — 1 (F — > 1) more rapidly than p+ approaches 
to zero. However, in this case, the term [[2]] causes a very 
strong repulsive effect since f{V) > and A ~ 1, while 
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the term [[3]] becomes negligible. This behavior can be 
analyzed in more detail as follows. Let us set V = 1 — S 
(S > 0). Then A 1 + (e'>'+ - 1)S + 0(6^), and one finds 
[[1]] ~ ^, [[2]] - S-\ and [[3]] ~ 6'/^ so that 

P+P+ ^ p+ 

p\ ~ eT - 1 ^ 

Therefore, when 5/ p+ 0, then the term [[2]] causes a 
very strong repulsive effect. 

In summary, the shell never reaches p+ = due to the 
strong repulsive force caused by the term [[1]] (when V 
is not so large), or by the term [[2]] (when V ^ 1). 

To make the statements below as concise as possible, 
let us first define the "core-region" . Consider the phase 
space r := {{p^, p^)\p^ > 0, — 1 < /j+ < 1}. It is under- 
stood that suitable topology is endowed on r.[13] At each 
time t, the region C is defined such as (i) [[2]] + [[3]] > 
and /?_)_ < on C and (ii) C is connected, and contains 
a neighborhood of (/?+ = 0,/)+ ~ ^1)- In other words, 
C consists of phase points that are "close enough" to the 
singularity-emerging point (p+ =0,p+ = —1). 

Definition 

We call the above-mentioned region C in F the core- 
region (at t), for brevity. 

We also say that the shell "enters (or, leaves) the 
core-region at t = t^," if at i = — the shell's 
phase point is outside (or, inside) the core-region, 
and at t = f*, the phase point of the shell is on 
the edge of the core- region; ([[2]] -|- [[3]] = and 
P+ < 0). 

It is clear (from the consideration of the order of mag- 
nitude of the terms [[1]], [[2]] and [[3]] when F ~ 1) that, 
if the shell would have ever reached the singularity, the 
phase point of the shell should have reached the point 
{p+ = 0,p+ = — 1) through the core-region. 

Now we can show 

Lemma 2 

Under the dynamical equation Eq.(35), once the 
shell enters the core-region at t = t*, it bounces 
back without reaching zero-radius. Indeed, p+{t) is 
bounded from below as 

P+{t) > P+{t*)\J 1 - P+{t*f . 
Proof: 

^_ .2 

Noting that p+ > — when the shell is 

in the core-region, it suffices to show the 
claim holds even for 



with arbitrary initial conditions p(0) = 
a{> 0) and p{0) = —b (where a > and 



1 > 6 > 0). However, its solution is given 
by p{t) = {{t - abf + W{1- 62))V2, so 
that p+{t) > p{t) > aVl - 6^ > 0. Thus 
the claim follows, q 

Theorem 1 

Under the dynamical equation (with or without Pq) 
Eq.(35), the shell never reaches p+ = irrespective 
of its initial conditions. 

Proof: 

It suffices only to consider the contract- 
ing phase. Then if the shell would have 
ever reached p+ = 0, the phase point of 
the shell should have reached the point 
{p+ = 0,p+ = —1) through the core- 
region. However, Lemma 2 indicates that 
this never happens, q 

Thus the line-like singularity never forms in this model 
irrespective of the value Pq, as far as the weak energy 
condition is satisfied outside the shell. 

The behavior in the expanding phase (and/or after the 
bounce) varies depending on the value Pq. When Pq>0, 
p+ > in the expanding phase, so that the shell contin- 
ues to be accelerated and expand, and its velocity asymp- 
totically approaches to the light velocity. On the other 
hand, the case Pq < 0, the term [[3]] is negative, thus 
it might be possible that p turns to negative when the 
shell velocity is not very large. Therefore, the oscillating 
behavior can arise, which might be interesting to investi- 
gate further. (However, even in this case, the shell never 
collapses to zero-radius due to Theorem 1.) Once the 
shell velocity becomes sufficiently large in the expand- 
ing phase, however, the shell continues to be accelerated 
and expands forever just as the case of Pq > 0. This is 
because, when the shell velocity becomes close to 1, by 
setting p+ = 1 — (5 ((5 > 0, (5 <C 1), the term [[2]] behaves 
as ~ and become dominant, while the terms [[1]] and 
[[2]] behave as ^ 5. 

Now let us illustrate what has been proved rigorously 
above by some numerical investigations. We set Pq = 
for simplicity. 

Figure 1 shows a typical behavior of the contracting 
shell (thinner curve) along with p^ (thicker curve). As 
an example, we have set 7+(7'+ — t^) = Y^(f'+ — i+ + 
100)'^. Initial conditions have been set p+(0) — 0.1 and 
p+(0) = —0.999 (i.e. very close to the fight velocity). 
The bouncing behavior discussed above is obviously seen 
in the figure. It is also seen that the turning point of the 
shell coincides with the peak of the rotational pressure. 

Figure 2 shows the behavior of the same model with 
the initial conditions /9+(0) = 0.1 and p+{Q) = —0.1 (i.e. 
relatively mild contraction). Again, evolutions of p+{t) 
(thinner curve) and p^j, (thicker curve) are indicated. Due 
to the mild initial conditions, the shell soon bounces 
back. Just as the previous case, the turning point of the 
shell coincides with the peak of the rotational pressure. 
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P, lOOr 




FIG. 1: Typical evolutions of the shell-radius (thinner curve) and the rotational pressure (thicker curve). We have set 
7+(r+ — t-f-) = ■^{r+ — t+ + 100)^. Initial conditions are p+(0) = 0.1 and p+(0) = —0.999. The vertical line indicates 100p+ 
and P0, while the horizontal line indicates t+. 




FIG. 2: The same model as FIG. 1 with initial conditions, p+(0) = 0.1 and p+(0) = —0.1. The vertical line indicates 10p+ (for 
the thinner curve) and (for the thicker curve), while the horizontal line indicates t+. 



Noting that [Kg^eJ = 0, Eq.(19) with Pz = implies 
e = Pij,, so that 



-7+ 



Ke = Kp^ - KPz = 



(A-1) 



np^, 



e-7+(A- 1) 



(37) 



This relation provides another confirmation of the valid- 
ity of the condition 7 > (sec Eq.(34)). 

Finally let us note that the relation Eq.(37) reveals a 
curious character of the model. We first consider the case 
Po = 0. In the case of the counter-rotating dust cylinder, 
there is a relation ^ = u^, where v is the velocity of 
the dust particle relative to an observer "standing still" 
on the shell [5]. Comparing this with Eq.(37) (Po = 0), 
the present model corresponds to the case of w = 1, so 



that the present model could be interpreted as a cylinder 
shell consisting of the counter-rotating massless particles, 
though the two models are not exactly the same. The sit- 
uation in which massless particles are confined to a com- 
pact region arises naturally in superconducting objects, 
so that the present model should have some relevance. 
On the other hand, when Pz < 0, it behaves more or less 
like a normal matter and it is consistent with the attrac- 
tive nature of the term [[3]]. On the other hand, when 
Pz > 0, Eq.(37) rather looks like a cosmic string. How- 
ever, the cylinder tends to expand in the present case, 
contrary to the case of cosmic strings. It may not appro- 
priate at present to pursue such comparisons too further. 
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IV. COMPARISON WITH THE RESULTS OF 
THE PRECEDING PAPER 

In view of our results above, we need to investi- 
gate whether the singularity-forming solution claimed in 
Ref. [6] and its correction Ref . [7] is a valid one. It turns 
out that the formal solution in question is not a relevant 
one and should be discarded. 

By comparing one by one our formulas with the corre- 
sponding ones shown in Refs.[6] and [7], it turns out that 
some formulas in Ref. [6] along with its correction Ref. [7] 
still need slight modifications (sec below and Appendix 
D for more details). However, we first derive the above 
conclusion faithfully based on the formulas in Refs.[7] and 
[6]. Then we also show that the same conclusion follows 
based on our formulas, too, even with the Pg-term. 

Firstly their functions /, g, h, I, Ro{T), Ro'{T), 
Ro"{T) and b{T — R) correspond to ours as 

f ■.= T^ , g := i?2 , h:=Z^ , I := $2 , 

i?o(r) ■.= p{t) , i?o'(T) ■.= p{t) , 

i?o"(r) :=p(t) , 6(T-i?) = -27(t-r) . (38) 

In particular, let us note that the sign of the function h is 
opposite to our function 7. Thus we should understand 

that & < and 6' < 0. They correspond to the physical 
conditions 7 > and 7' > (see arguments after Eq.(24) 
and Eq.(25)). 

The detailed comparison is shown in Appendix D. 
Here it suffices to pay attention only to their dynami- 
cal equation shown in Ref. [7] (which is the correction of 
the original Eq.(30) of Ref. [6]): 



(i?o' - A) 



^.2 [ A y(Co)(l-j?oO 
i?o -(l-i?o )|^+ 2(A-1) 

(39) 

where Rq and its derivatives correspond to our /9_|_ {t) and 
its i-derivatives (see Eq.(38)), and 6(^0) corresponds to 

our -27+; A := [i?o'^ + e-\l - i?o'')] ^'^ ■ 

If the terms are rearranged to look as similar to Eq.(35) 
(with Pq = 0) as possible, their formula Eq.(39) reads 



i?n" = A 



I--R0 



/2 



+ 



i?0 



2(A-1) 



(1 - i?o' )(1 - i?o')(A - Ro') . (40) 



(Compared with Eq.(35) (Pq = 0), the factors (1 — Rq'^) 
and (A — _Ro') of the second term should be replaced by 
i?o' and ((A + l)Po' + l), respectively.) 

Based on Eq.(39), Ref. [7] (just as Ref. [6]) tries to 
search for special solutions for Rq satisfying 



Ro" = 

a ■= Ro' 



1 - Ro" 
Ro 

A {-b'){l-Ro') 
Ro 2(A-1) 



(A - Ro'] 



(41) 



, (42) 



where /? is assumed to be an arbitrary constant. 

As is claimed in Ref. [6], Eq.(41) itself is formally sat- 
isfied by 



i?o'(r) = -(l-e-2/3(^-^°))i/2 



(43) 



or integrating once more, 



Ro{T) = {R,-{T-To)} 

+ i|(l _ g-2/3(T-To))l/2 

-ln[l + (l-e-2'5(^-^«))i 



))V2]} .(44) 



Here two integral constants Tq and i?* have appeared 
corresponding to the fact that the differential equation 
Eq.(41) is of the 2nd order. The constant Tq indicates 
the time when Rq vanishes, while i?» (along with Tq) 
determines the time T* when Rq vanishes (T* < Tq). 
These constants To and i?,* are (indirectly) determined 
by the initial conditions for solving Eq.(41). 

Equation (44) would imply a line-like singularity for- 
mation, reaching i?o = at some finite time T*, so that 
its implication would be significant. We show below, 
however, that Eq.(4-4) is not the solution of the original 
differential equation Eq.(40), irrespective of the parame- 
ters To and R^, , i.e. independent of the initial conditions. 

To see this conclusion, let us first define B{T) to be 
-1 times the R.H.S. of Eq.(42), 



B{T) := i-Ro') 



A 
R^ 



i-b'){l-Ro') 
2(A-1) 



{A-Ro') 



(45) 

What we show below is that B{T) varies (actually di- 
verges) along the collapsing solution Eq.(44) so that the 
starting assumption /3=constant (Eq.(42)) is incompati- 
ble with the dynamical equation Eq.(41). Causing such 
self-inconsistency, then, Eq. (44) is not the solution of the 
original dynamical equation, Eq.(40). 
Now we prove 

Theorem 2 

The quantity B{T) diverges, and does not remain 
constant along the collapsing solution Eq.(44)- 

Proof: 

Let us look at the definition of B{T) 
(Eq.(45)). Both of the terms in the brack- 
ets [ ] are positive, since — 1 < Rq' < 0, 
6' < and A > 1. Thus B{T) > and 
accordingly, Eq.(42) implies that (3 < 0. 
We then get the estimation 

B{T) > i-Ro')^ > t^ol , 

SO that 

lim BiT) > lim 

T^T* T^T* Rq 
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where T* is the shell-collapsing time, 
limx^T" RoiT) = 0. 

In the case lim-r^T* Ro'{T) ^ 0, the 
above estimation imphes that B{T) — > oo 

as T ^ T*. 

In the case Iiuit^t' Ro'{T) = 0, the 
estimation further goes as 

lim B(T) > hm ^ 

T^T» T^T* Rq 

i-Ih") 

where de L'Hopital's theorem has been 
applied in the last step. With the help 
of Eq.(41), then, we get 



lim B(T) > (-13) lim - — ^MlL 



= lim 



> ^(-P) lim — 5- 



oo 



Thus we get lim^^T* B(T) = oo. |— | 

The situation is the same for our correct equation, 
Eq.(35) (with Pq > 0, = or < 0), too. Just as Eqs.(41) 
and (42), wc can ask whether the collapsing solution is 
allowed for the equation 



P+ 
r A 

(3: = P+[— + 



(46) 



7P+ 



A- 1 



— {(A + l)p+ + l} 
(l-p+2)V2eT+Pol , (47) 



p+ (A-l)(l + p+) 
A 



with [3 being a constant; A := 
[(-p+)2+e2''+{l-(-p+)2}]'/'. Equation (46) it- 
self has a formal solution, which is just the same as 
Eqs.(43) and (44). Though more care is required, 
essentially the same argument as above follows and it 
turns out that the formal solution just like Eq.(44) is not 
the solution of the original dynamical equation Eq.(35). 
To show this fact, we first introduce some symbols. 

Definition 



(a) Let B be -1 times the R.H.S. of Eq.(47); 

]3{t) := y X 



Lp+ (A-i)(i-y)' 

where V := -p+ (0 < < 1). 



(48) 



(b) Let D and E be the second and the third terms in 
the brackets [ ] in Eq.(48), respectively; 



m ■■ = 

m : = 



(A-m-v) 

A 



[(A + 1)F-1] 



1 



(I_y2^1/2e7+p(j . 



(c) Let /(V) := (A + 1)V - 1. The sign of /(V) deter- 

mines the sign of D{t). As is already discussed in 
the previous section (after Lemma 1), the function 
f{V) is monotonic, increasing, continuous on [0. 1], 
with /(O) = -1 and /(I) = 1, so that f{V) has only 
one zero, say V* , in (0,1). When D is looked as 
a function of V, then, D < for F e (0, V*] and 
D > for 1/ e {V*,l). 

(d) Let t* be the shell-collapsing time (which is assumed 

to exist in the present context), i.e. t* is the time 
such that limt^j. P+it) = 0- (It does not matter 
whether t* is finite or infinite in the argument be- 
low.) 

We now prove a theorem {Theorem 3 below) corre- 
sponding to Theorem 2 step by step. 
We first note 

Lemma 3 

\E\/iD^Q asV 

Proof: 

When V is close to 1, setting V = 1 — 5 
{5 > 0), then A = l + (e'>'+ - 1)6 + 0(^2), 
and one finds 1/\D\ = 0(5^) and 1/\E\ = 
0((5i/2), SO that \E\/D = 0(^3/2). Then 
the claim follows, q 

Now let us confirm 

Lemma 4 

It follows that limj^f V{t) < 1 . 

Proof: 

Studying the solution of Eq.(46) (cf. 
Eq.(43)), we see V < 1 (i.e. p+{t) > -1). 
Suppose limt^j. V{t) = 1. Since f{V) 
is continuous on [0, 1], then, it means 
limt^t* fiV{t)) = 1, so that D{t) oo 
as well as ^ oo as t ^ (note A > 
1). Considering the limiting behavior of 
B{t) = V{t)[^+Dit){l + Eit)/i){t))]BS 

t ^ t*, Lemma 3 implies B{t) oo as 
t — > t*. Since /3 (= —B{t)) is a constant 
by assumption, this should not happen. 
Thus limt^t. V{t) 7^ 1. □ 

Corollary of Lemma 4 

The quantities D{t) and E{t) are bounded, \D{t)\ < 
oo and \E{t)\ < oo. Furthermore, limt_>(. \D{t)\ < 
oo and limt^j. < oo. 
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Proof: 

We note that D and E become unbounded 
only when V 1 and that D and E 
are continuous as functions of V . Now, 
the solution of Eq.(46) (cf. Eq.(43)) im- 
plies y < 1. Furthermore Lemma 4 says 
limt^t. V{t) < 1. Then the claim follows. 

□ 

Then we can show 
Lemma 5 

The constant parameter $ should be set negative, 
Proof: 

We note that limj^t. < oo and 

limt^j. \E{t)\ < oo due to Corollary 
of Lemma 4- On the other hand, 
limt^f ^ = 00. Then one can find ap- 
propriate ti and t2 {ti < t2 < t*) such 
that A > \D{t)\ + \E{t)\ on h ■= ih, t^). 
Thus B{t) > on /i. Since j3 (= -B{t)) is 
a constant by assumption, this means that 
(3 should be set negative from the outset. 

□ 

Finally we prove 
Theorem 3 

The quantity B{T) diverges, and does not remain 
constant along the collapsing solution (correspond- 
ing to Eq.(44)). 

Proof: 

Let us set C{t) := B{t) - {-p+)[D{t) + 
m]. We note C{t) = (> 0) 

due to Eq.(48). Then it follows C{t) > 
^^^+i£)i. Noting /3 < {Lemma 5), then, 
exactly the same argument as the proof of 
Theorem 2 can be applied, and one con- 
cludes limt_»t. C{t) = 00. On the other 
hand, it follows 

lim C(t) < lim \B(t)\ 

+ ^lim \{-p+)Dit)\ + lim \{-p+)m\ . 

Due to Corollary of Lemma 4, then, it fol- 
lows limt_>(. \B{t)\ = 00. □ 

We conclude that there is no relevant solution indicat- 
ing the singularity formation in the present model. 



SUMMARY 



We have started with constructing a general framework 
for analyzing a cylindrical spacetime and a shell in it, 
which might be useful for future investigations. Based on 
the framework, we have investigated a cylindrical shell- 
collapse model which accompanies the out-going radia- 
tion of gravitational waves and massless particles. This 
model had been introduced by Pereira and Wang [6, 7], 
but its proper analysis had been awaited. This model 
could be interpreted as a thin shell filled with radiation. 
In connection with the cosmic- censorship hypothesis, we 
are mostly interested in the collapse of "normal" mat- 
ter, so that the weak energy condition has been assumed 
outside the shell. 

It has been proved that, as far as the weak energy 
condition is satisfied outside the shell, the shell bounces 
due to the rotational-pressure effect. After the bounce, 
it continues to expand without re-contraction when the 
pressure of the shell in the z-direction pz satisfies Pz > 0, 
while in the case Pz < 0, the behavior after the bounce 
can be more complicated. However in either case, the 
shell never reaches the zero-radius configuration and it 
escapes from the line-like singularity formation. 

Just as the case of a shell filled with counter- 
rotating dust particles considered by Apocatopolas and 
Thorne [5], the present case also shows a bouncing due to 
the effect of rotational pressure, while the eternal expan- 
sion after the bounce (when pz > 0) is a unique feature 
of the present model. We have also performed numerical 
investigations which reveal explicit behaviors of the shell. 
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APPENDIX A: SUMMARY OF 
FUNDAMENTAL GEOMETRICAL QUANTITIES 
FOR CYLINDRICAL SPACETIMES 



In this paper, we have studied the dynamics of a cylin- 
drical shell in the spacetime cylindrical symmetry. 



Based on the metric Eq.(l), it is convenient to intro- 
duce the 1-forms as Eq.(2). Taking the exterior derivative 
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of ^ s, we get 

TR ' TR 



(Al) 



Here " ° " indicates the derivative w.r.t. the proper-time 
r of an observer on the shell, and X := dt/dr; they are 
evaluated on the shell in consideration. 



Taking the exterior-derivatives of w'^b's, we get 



Here " • " indicates the partial derivative w.r.t. t, while 
" ' " means the same w.r.t. r. With the help of the first 
Cartan structure-equation dO'^ +uj"^b =0, we thus 
get 

, ,0 — fl" ^ fll , — ^ (fi 

, ,0 _ * fl3 , ,1 _ fl2 



(A2) 



On account of the relation, w^b = ^^bc we also 
get [14] 



i 22 — 



1 o o , 

iz-pz') , 



TZ XTZ 

$ loo 



r"33 = — , 



r^22 = — 



Z' 



RZ 



r^33 = — 



others = . 



(A3) 



d(jj°i = 

diO°2 = 

dw°s = 

dLV^2 = 

du'^'i = - 



'r\ /T'Y 



TR I \T 

1 / 



R 



9° A 



A6»3 



7'\' 



1 ^*'UoA«3 1 m'^..^^ 



V R 



duj'^3 



= 



(A4) 



Thus, by the second Cartan structure-equation VI^b = 
du^B + i^^c A LO^ B, we find 




TZ \T R"^ TZ \ ] RZ \T TR RZ 
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Using the relation Q^b = ^R-scd ^'^ ^ ^® obtain the Riemann curvature w.r.t. the frame {e^}, 



1 T'Z' 




R'°212 = 



RZ 



T 



1 RZ' 

tr'rz 

1 T'Z 



J_ lZ^\ 

'tz \ r] ^ trtz 



313 



1 



Thus, the Ricci curvature in the frame {e^} becomes 



■00 



■01 



Rii = 



R 



22 




1 R 

TRR 



Z^ 



- 1 

R I T 



1 



1 



^1 




1 Z' 






RZ ' 




^ R^Y 







R33 = 

others = 

Finally, the scalar curvature becomes, independently of the frame. 



J-l 




' 1 






R^ ' 




" ^¥ ' 


T" 





R 



2 l/^\ l/| 



1 /T' 



^1 




1 T' 






TR ' 












1 



(A6) 



(A7) 



(A8) 



APPENDIX B: NON-COORDINATE BASES 



{ca} (see Sec. II A), 



Throughout this paper, we always consider a cylindri- 
cal spacetime defined by the metric Eq.(l). The coordi- 
nates {x^} :— {t, r, z, (/>} define the coordinate bases 
and their dual 1-forms, 



Sfj, '■ — dfj, , 



r := dx^" 



(Bl) 



(Here the Greek letters /x, i^, • • • indicate the indices w.r.t. 
the coordinates {t, r, z, <p}.) 

On the other hand, a cylindrical shell embedded in the 
spacetime naturally defines a local ortho-normal frame 



^ p, Ix, 0, 

/ {trz<i>) 

ir = dr = (X, P, 0, 0)(t,,^) , (B2) 

= (0, 0, 1/Z, 0)(,,,^) , 

= (0, 0, 0, l/<I')(t„^) , (B3) 

where the symbol e implies the normalized frame- vector. 
(The first few Greek letters a, /3, 7, • • • are used for the 
indices w.r.t. the above ortho-normal frame, and they 
take the values n, t, z or 0.) A set of 1-forms {9"} dual 
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to {Ca} is also introduced. 

As the third set of frames, the metric given in Eq.(l) 
defines a natural set of 1-forms {0^}, 



e° = T{t,r)dt , = R{t,r)dr , 
192 = Z{t,r)dz , 9^ = <!>{t,r)d(j) , 



(B4) 



along with their dual ortho-normal bases {ca}- (The 
capital Latin letters A, B, C, ■ ■ ■ indicate the indices 
w.r.t. the frame determined by Eq.(B4), and they take 
the values — 3.) 

Now we summarize the relations among {en}, {ca} and 
{e^}, and those among {6"^}, {O"} and {0^}. 

Firstly, Eq.(B3) gives the relation between {e^} and 
{ia}, and those between {9'^} and {§"} 



(B5) 



Here E^a and E^^ defined below are mutually inverse as 
matrices: 



T-^Rp° X 

XTR-^ p° 





/ -TRp° XTR \ 

-i?2p° 

Z 

V $y 



(B6) 



Secondly, Eq.(B4) gives the relation between {e^} and 
{e^}, and those between {9^^} and {9^} 

9^ = E^i, r , 9" = E^a9^ , 

BA = E^A , en = eA E^^ . (B7) 

Here E''*^ and Ei^a defined below are mutually inverse 
as matrices: 

E^n = diag{T, R, Z, $) , 

E^'a = diag{l/T, 1/R, 1/Z, 1/$) . (B8) 



Finally, Eqs.(B5) and (B7) give the relation between 
{ca} and {ia}, and those between {9^} and {9"} 



e« = caE^ 



= E"a9^ 



(B9) 



ReveE^a — '^'^ 
are mutually inverse as matrices: 



E^^E^^a and E"" a ■= E'^^EI'a given below 



E' 



E"a = 



/ Rp° XT 

XT Rp° 

10 

V 1 

-Rp° XT 0' 

XT -Rp° 

10 

1 
Now noting the relation [12], 



(BIO) 



-^7 ^ 



along with Eq.(B9), we get 

= E^„E''0E'^cT%+E<'o.E^BdnE''0.{Bll) 

Equation (Bll) is the basis for calculating the extrinsic 
curvature. 



APPENDIX C: DERIVATION OF THE 
FORMULA FOR /sTe^e, 

Here derivations of the formula for the extrinsic cur- 
vatures (Eqs.(13) and (16)) are purposefully shown in 
detail. This is to ensure anyone interested in this topic 
reproducing all the results presented in this paper with- 
out any difficulty and to motivate him to go ahead the 
point where this paper ends. 



In order to derive Eq.(13), we note that the first term on the R.H.S. in Eq.(12) becomes, with the help of Eqs.(A3) 
and (BIO), 

jpA TTiB zpn _T^C TTil TTiO rpn j^O /' zt'I \2 rpn t^O 

— r-C/ T-C/ Ci- AB — ~^ T-C/ r-C/ O-L lo ~ r) 0^- n 



= {RP)- 



( XT ,^ R o 

< (InT^)' + p 

\2R^ ' 2XT 



{In Hy - p{\nR^y 



(CI) 



Similarly, with the help of Eqs.(B6) and (BIO), the second term on the R.H.S. in Eq.(12) is calculated as 
—E^T-E"' Bdij,E^ T = —E^o {E^T-dtE^T + E^T-drE^T-) ~ E"'i {E^T-dfE^T + E^T-drE^T^ 
= Rp{XT)° - XT{Rp)° . 



The above expression can be further modified with straightforwardly; 

the help of the following formulas that can be derived o y d2 2 

^ - ^ p p-^{l^Tr + ^,P (lni?^)° 



XT^ 2 ^ ^ 2XT2 

X' = -lj^(lnT7 + ^p (lni?2)' . (C2) 



16 



After some calculations we then get 

1 



2\o 



-P (InT^)'- 
2 2XT 

-. {I) + {II) + {III) + {IV) . 



The last two terms together yield 
(7/7) + {IV) 



XTR o2 



2 /d2\o 



XTR a„(lnT77?^) - R^ P {\tiT''/R 

2XT 

=■■ {Vi) + {V2) + (F3) + (Vi) , 



R 



XTR 9„(lni?2) _ i?2 p (ini?2)<: 



where the formulas in Eq.(lO) have been used. Noting 
that (77) + (14) = -2§tP ° (ln-R^)°, thus, the calcula- 
tions go as 

Ke^e. = (7) + (77) + (777) + (7F) 

= {I) + {iII) + (y,)} + {Vs) + {{V) + {V,)} , 

yielding Eq.(13). 

Now we derive Eq.(16). To start with, let us note the 
following formula derived straightforwardly from Eq.(7): 



-(InT^ + 



p^(ln7?^) +X^R' pp 



Taking the i-derivative of the relation P ° = Xp (Eq.(7)), 
then, we obtain. 



1 00 1 o „ i?2 o3 



We should return to the basic equation Eq.(12) along 
with the calculations Eqs.(Cl) and (C2). Noting the re- 
lations 

{XT)° = X{{TX) + p{TX)'} 
XT 

= —{X^T^-l){\n{R^/T^)y+X^TR^pp , 

o o o 00 

{Rp)° = Rp+Rp 

= X^T^Rp+^X^T'^R p {lji{R'^/T^))° , 



(C3) 



we finally reach Eq.(16). 



APPENDIX D: DETAILED COMPARISON 
WITH THE FORMULAS IN THE PRECEDING 
PAPER 

1. Extrinsic curvature 

Here we point out the discrepancy of the formulas re- 
ported in Ref. [6] and its correction Ref.[7] with our for- 
mulas. 

Let us note once again that their functions /, g, h, I, 
Ro{T), Ro'{T), Ro"{T) and b{T - R) correspond to ours 

as 

f:=T\g:=R\ h -.^ Z\ I := $2 , 

Ro{T) ■.= p{t) , i?o'(T) ■.= p{t) , 

i?o"(T) := p{t) , b{T -R) = -2j{t - r) . (Dl) 

They show the expression for the (Tr)-component of 
the extrinsic curvature (see the first formula in Eq.(lO) 
of Ref. [6]), which we now denote KTr^\ as 



1 



ifg) 



1/2 



2(/-5i?o'(T)')3/2 

f 



x<-/iii + f/f -2^)i?o'(r) 



(2^-^) (i?o'(T))2 + ^i?o'(T)'-2i?o"(r)| . (D2) 



r 



It is easily checked that their result is translated into our 
notation as 



X'TRp+^ilnR') +^ilnTy 



-^X'^TR p {\n{T^/R^)y 



(D3) 



where the r-derivative in the last term is for making the 
expression concise, and can be expressed in terms of " • " 



and " ' " if needed. Comparing Eq.(D3) with Eq.(16), 
we realize, 

7^™ = -K^^^^ + ^{(InT^)' + ^ P {InR^i } 
= -Tfe^e. + ^^a„ln7?2 + _^a„ln(TV7i2) 



2 

X^T^ 



Rp{\n{T^/R^))° 



(D4) 
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Here we have made use of Eq.(lO) to reach the last line 
in Eq.(D4); 

= -B? P{HT'^/R^))° + XTR dn HT^/R^) 
+^5„lni?2 . 

On the other hand, their expressions for the (zz)- 
component and the (0(^)-coniponent of the extrinsic cur- 
vature, Kz^^^ and K^j^^^^ , (see the second and the third 
formulas in Eq.(lO) of Ref.[6]) are related to our Ke^e^ 
and Kej,eA as 



1/2 



2 \f-gRo'{T)\ 



(D5) 



K 



(PW) _ _1 



fg 



1/2 



(D6) 



It turns out that ifif and k'^'^^ are totally consis- 
tent with our Ke^e^ and Ke^e^, taking into account that 
the former is based on the coordinate frame {9^}/i=t.r,z,0) 
while the latter is based on the ortho-normal frame 
{e-a}a=n,T,z,4,- [Ou the other hand, Krr^^ is the (rr)- 
component (and not (tt)-componcnt), and dr = e,- is a 
member of the ortho-normal frame {ea}a=n,T,z,4>- There- 
fore K^T^^ is what should correspond to our ife^e^ (uP 
to sign) without any further conversion factor.] 
Looking at Eqs.(D4)-(D6), we realize 

(0) Their definition for the extrinsic curvature and ours 
differ by the overall sign, which is just the matter 
of definition and is not essential. 



(1) The last three terms in the last line on the R.H.S. in 

Eq.(D4) are the essential discrepancies in the (rr)- 
component of the extrinsic curvature. 

(2) There is no essential discrepancy regarding the (zz)- 

component and the (00)-component of the extrinsic 
curvature. 

Their expression Eq.(D2) can be modified so as to co- 
incide with our result Kg^g^; 

j^{PW,new) _ 



1 



ifg) 



1/2 



{f-gRo'{T)y/^ 



+ ^Ro'{Tf -2Ro"{T)[. 



(D7) 



Comparing Eq.(D7) with Eq.(D2), we realize that the 
term in Eq.(D2) should be removed and that the 

numerical factor "2" in front of the factor ^ should be 

9 

corrected to "1". 



2. Formulas for the pressures pz and 

The expression for ^ (corresponding to our ^) in 
Ref. [6] has been corrected in Ref. [7] , and the new expres- 
sion coincides with our Eq.(27). However, the expressions 

for and r'^it) (corresponding to and p-{t-), 

respectively) in Eq.(28) of Ref. [6] are not mentioned in 
Ref. [7]; they actually should be corrected as our Eqs.(28) 
and (30), respectively. The expression for has been 
corrected in Ref. [7] as 



J 



e«'(«o)/2 

i^Pz = — — ^tsTtt: <1 (1 - A) 



A(l - Ro'y/^ 



---(1 — Rq' ) — i?o" 



^^(l-i?o')(l-i?o'')(i?o'-A) 



(D8) 



r 



However, it still needs modification (sec our Eq.(32)), 
which is equivalent to replace the last two factors (1 — 

Ra'^){Ra'-A) in Eq.(D8) with -i?o'[(A-hl)i?o' + l]- The 
new expression of influences the dynamical equation 
for the shell, since the latter is obtained by setting p^ = 



(They only consider the case Pz = 0). The rest is 
discussed in Sec. IV. 

The expressions p^ in Eq.(29) of Ref. [6] should also be 
modified as our Eq. (33). 
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